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Observation of Kosterlitz-Thouless spin correlations in the colossally magnetoresistive 

layered manganite Lai 2Sri 8Mn207 
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The spin correlations of the bilayer manganite Lai.2Sri.gMn207 have been studied using neutron 
scattering. On cooling within the paramagnetic state, we observe purely two-dimensional behavior 
with a crossover to three-dimensional scaling close to the ferromagnetic transition. Below Tc, an 
effective finite size behavior is observed. The quantitative agreement of these observations with the 
conventional quasi two-dimensional Kosterlitz-Thouless model indicates that the phase transition is 
driven by the growth of magnetic correlations, which are only weakly coupled to polarons above Tc- 
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It is now generally accepted that colossal magnetore- 
sistance (CMR) in doped manganese oxides involves a 
strong coupling between spin, charge, and lattice de- 
grees of freedom ||,|^ close to the ferromagnetic tran- 
sition. However, the nature of the transition and the 
relative importance of the different interactions in con- 
trolling the magnetotransport above the critical temper- 
ature Tc have not been clearly established. A number 
of studies have indicated that the magnetic phase tran- 
sition may be unconventional in CMR compounds, with 
the observation of an anomalous spin diffusion compo- 
nent below Tc and a non-divergent correlation length . 
This has been attributed to the development of magnetic 
polarons above Tc and their possible persistence below 
Tc, although it can be difficult to distinguish them from 
standard critical fluctuations. In other studies, conven- 
tional critical scaling of bulk properties has been observed 
Q. It is important, therefore, to study the nature of the 
phase transition in CMR compounds, and determine if 
unconventional magnetic correlations are essential to the 
mechanism of CMR. 

Naturally layered manganites are derived from the 
perovskite structure of the three-dimensional (3D) com- 
pounds by the addition of non-magnetic blocking layers 
|^,||. The bilayer compounds La2_2a;Sri+2a:Mn207, in 
which X represents the hole concentration on the Mn02 
planes, have been extensively studied in recent years be- 
cause of the insights they provide into the mechanisms 
of CMR @-|l|]. The motivation of the present study 
is to utilize the low-dimensionality of these compounds 
to perform a detailed investigation of the spin correla- 
tions close to Tc using neutron scattering. The reduced 
dimensionality of the spin fluctuations extends the tem- 
perature region over which critical fluctuations may be 
studied, and makes them easier to distinguish from other 
dynamic processes in the sample . This has allowed us 
to compare the temperature evolution of magnetic corre- 



lations with other two-dimensional (2D) systems exhibit- 
ing similar magnetic ordering. 

In previous work on the 40% hole-doped bilayer sys- 
tem we observed strong 2D in-plane ferromagnetic fluc- 
tuations above Tc, with evidence of competing ferromag- 
netic and antiferromagnetic interactions perpendicular to 
the planes . The in-plane correlation length ^ was mea- 
sured in scans through the 2D rods of magnetic scatter- 
ing chosen to optimize the energy integration. However, 
these measurements are necessarily performed away from 
the wavevector corresponding to 3D magnetic ordering, 
and so were not sensitive to a possible crossover to 3D 
correlations. 

We have now investigated the correlations close to the 
3D ordering wavevector as a function of temperature, 
both above and below Tc ■ We show that the spin corre- 
lations are quantitatively consistent with a quasi-2D XY 
model, exhibiting Kosterlitz-Thouless correlations above 
Tc with a crossover to 3D correlations at a correlation 
length consistent with the known in-plane and inter- 
bilayer exchange constants. This agreement with other 
2D systems indicates that the magnetic correlations de- 
velop conventionally, and that the transition is a true 
second-order phase transition, although there is evidence 
that it is smeared by weak inhomogeneity. This is in 
contrast to a magnetic polaron model, in which the mag- 
netic correlations are strongly bound to charge degrees 
of freedom. Recent x-ray and neutron scattering results 
provide evidence of charge localization and the develop- 
ment of charge correlations above Tc Nevertheless, 
they do not have a strong influence on the magnetic cor- 
relations above Tc- Instead, the incipient charge ordering 
is preempted by the onset of ferromagnetic ordering, and 
the charge correlations collapse at Tc- 

The neutron scattering experiments were performed at 
the NIST Center for Neutron Research, using the same 
single crystal of the 40% hole-doped bilayer manganite 
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FIG. 1. Magnetic diffuse scattering of Lai,2Sri,8Mn207 
measured along Q — [h02]. The solid lines are the results of 
the lineshape analysis, including a Gaussian Tc-distribution 
and the instrumental resolution, as explained in the text. The 
dotted line for T = 113.5 K denotes the best least squares fit 
to a Lorentzian profile convoluted with the instrumental reso- 
lution function. The inset shows the temperature dependence 
of the (002) Bragg peak from Ref. Q and the result of a fit to 
a standard power-law behavior convoluted with a Gaussian 
Tc-distribution. 

Lai.2Sri.8Mn207 (lattice parameters a = 3.862 A and 
c = 20.032 A at 125 K) that has already been studied in 
various other experiments [p|-pl|. This compound shows 
a transition to long-range ferromagnetic order with the 
Mn spins aligned entirely within the a-b plane at Tc ~ 
113 K, coincident with the metal-insulator transition ob- 
served in resistivity measurements. Details on the prepa- 
ration and characterization of the sample are given by 
Mitchell et al. 

The magnetic diffuse scattering was measured on the 
BT2 triple-axis spectrometer operating in two-axis mode, 
i.e., without analyzing the energy of the scattered neu- 
trons. Scans were taken around the (002) Bragg position 
where the nuclear scattering contribution is extremely 
weak, with a fixed incident energy 13.7 meV and horizon- 
tal coUimations of 60'-20'-20', full width at half maximum 
(FWHM). Pyrolitic graphite was used both as monochro- 
mator and filter against higher order contamination. In 
this two-axis mode, the diffuse scattering is, in the qua- 
sistatic approximation, proportional to the wavevector- 
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FIG. 2. Temperature dependence of the correlation length 
^ obtained from a Lorentzian profile analysis of the scans 
along Q = \Ml2] (circles) and Q = [l+h, 0, 1.83] (open squares, 
from Ref. g ). The solid line denotes a fit to the Q2DXY 
model as explained in the text. 

dependent susceptibility XT(q), where q = Q — r, Q is 
the momentum transfer and r denotes a reciprocal lattice 
vector of the magnetic structure. 

Figure |l] shows the magnetic diffuse scattering ob- 
served in scans along Q — [/i02]. The susceptibility is 
well described by the usual Lorentzian profile except for 
temperatures close to Tc- In this region, however, there 
is also a strong tail in the order parameter that was at- 
tributed to an inhomogeneous broading of the transition 
common to disordered systems like the doped mangan- 
ites . Because of the strong dependence of Tc on hole 
doping in the bilayer manganites [ p^ , such sample inho- 
mogeneities would manifest themselves as a distribution 
of Tc ■ We have therefore reanalyzed the scaling behav- 
ior of the order parameter, allowing for a Gaussian dis- 
tribution of Tc's. As the inset to Fig. Q shows, such a 
distribution is in excellent agreement with the observa- 
tions. The newly obtained values of the order parameter 
are (i = 0.14(1) and Tc = 113.2(2) K, compared to /3 = 
0.13(1) and Tc = 111.7(2) K reported in Ref. §. The 
derived standard deviation gtc = 1.6 K of the Gaussian 
T(7-distribution corresponds to less than 0.4% variation 
in the total hole-doping |lj] . This shows that even a very 
small sample inhomogeneity strongly affects the measure- 
ments and has to be included for the analysis of critical 
quantities. 

An analysis of the wavevector dependent susceptibility 
including the Tc-distribution requires a detailed knowl- 
edge of the phase transition, i.e., a model for the tempera- 
ture dependence of the correlation length ^ and the static 
susceptibility xt(0), which will be derived below. For 
the following discussion, we use the correlation lengths 
obtained from a Lorentzian lineshape analysis, which are 
shown in Fig. For these to be valid, the scans need to 
integrate over all fluctuation frequencies. In 2D systems, 
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an optimal energy integration is achieved for a scatter- 
ing geometry that has the scattered wavevector parallel 
to the c-axis [|l6|, but this is not kinematically possible 
in scans along Q = [h02]. However, the agreement with 
previous measurements performed in the optimal config- 
uration at Q = [1 + /i,0, 1.833] (Rcf. g) vaUdates the 
quasistatic approximation in the present experiments. 

In the layered manganites, the separation of the Mn02 
bilayers by insulating (La,Sr)0 layers leads to a strong 
anisotropy, both in transport properties and magnetic 
correlations. In a simple nearest-neighbor exchange 
Hamiltonian for the magnetic interactions, one therefore 
expects large differences between the various exchange 
constants: Ji between spins within the same plane, J2 be- 
tween spins in different layers within a bilayer, and J3 be- 
tween spins in different bilayers. This quasi-2D behavior 
has been verified by spin-wave measurements |pl] , p^ , p^ 
which yield J1/J3 ~ 150, similar to the anisotropy ob- 
served in the resistivity ||l5| . 

To our knowledge, the critical behavior of quasi-2D bi- 
layer systems has not yet been investigated thoroughly, 
but there exist detailed theoretical and experimental 
studies of single-layer quasi-2D XY magnets (Q2DXY), 
in which there is a strong exchange anisotropy J/ J', 
where J and J' are the in-plane and interlayer exchange 
constants, respectively []l7|-pOt. In such systems, both 
above and below Tc, spin fluctuations on length scales 
that are small compared to the characteristic length Lcs 
« \/{J/J') (in units of the nearest neighbor Mn-Mn 
distance a) are not affected by the interlayer coupling 
and are therefore purely 2D. When approaching Tc from 
above, one expects the correlation length and the static 
susceptibility XT(q = 0) to increase according to the 
Kosterlitz-Thouless expressions 



400 



and 



e/a = 6exp[5(r/TKT-l)-'/'] 



Xt(0) = Cexp[S(r/TKT - 1)"'/'], 



(1) 



(2) 



where Tkt is the topological ordering temperature of 
the 2D XY model |2^. Once the interlayer interac- 
tion becomes important, i.e., when the correlation length 
reaches the order of Leff, a crossover to 3D scaling is ex- 
pected. Renormalization group theory estimates a 3D 
ordering temperature Tc that is related to Tkt by |18[ 



(3) 



Below Tc, the correlation length decreases rapidly to Loff, 
where spin fluctuations are again unaffected by the in- 
terlayer coupling and the correlation length remains con- 
stant. From this point on, although the magnetization 
is 3D, the fluctuations arc 2D and the system can be 
modeled by a 2D system of effective size Loff |19|. This 
behavior has been observed, for example, in Rb2CrCl4 
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FIG. 3. Temperature dependence of the static susceptibil- 
ity xt{Q) (circles) and quasistatic polaronic scattering (di- 
amonds) from Ref. [^. The solid line denotes a fit to the 
Q2DXY model as explained in the text. The dashed line is a 
guide to the eye. 



As is shown in Fig. ^, our observations are in remark- 
able agreement with the above predictions of the Q2DXY 
model. From a least squares fit of Eq. (|l]) to the observed 
correlation length above 120 K, we obtain ^0 = 0.3(1), 
Tkt — 64(5) K, and b = 2.1(2), in excellent agreement 
with the theoretical value of 6 ~ 1.9 ||l^,|2^. For the 
static susceptibility we obtain B — 3.9(4). very close to 
the theoretical value B = b * {2 — j]), where the critical 
exponent rj < 1/4. Using the fitted values for b and Tkt 
and J/ J' = 150, we obtain Tc ~ 109 K, consistent with 
the value Tc = 113.2 K derived from the order parame- 
ter. 

The deviation from the purely 2D behavior above Tc 
occurs at a correlation length ^2D ~ 12 A, which also 
corresponds to the effective finite size observed below Tc- 
This is in rough agreement with the value of Tcff ~ 40 
A predicted from the measured in-plane and interbilayer 
exchange constants. In the region around Tc, we approx- 
imate the behavior of both the correlation length and the 
static susceptibility with an abrupt crossover to standard 
3D scaling, although a more complete description would 
have to include a crossover scaling function |^ . 

This model for the correlation length and the static 
susceptibility (shown as solid lines in Figs. || and ^ ) can 
now be used for an analysis of the peak lineshapes that 
includes the Tc-distribution. Figure |^ shows that this 
model gives a good description of the lineshape even close 
to Tc- It is important to note that the model includes 
a divergence at Tc of both ^ and xt(0), which is not 
apparent in the observed data. From a simple Lorcntzian 
lineshape analysis, one would conclude erroneously that 
the correlation length does not diverge at Tc- 

Although the inhomogenous broadening prevents a re- 
liable determination of the critical scaling within the 3D 
regime, it is evident from the above analysis that the 
magnetic correlations are quantitatively consistent with 
a conventional model of quasi-2D behavior. To underline 
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the significance of this conclusion to our understanding 
of CMR, we compare our results with the recent obser- 
vation of charge correlations in the same compound 
In the paramagnetic phase, we have observed a growth 
of diffuse x-ray and neutron scattering, arising from the 
strain field produced by quasistatic polarons. Further- 
more, broad incommensurate peaks modulating this dif- 
fuse scattering show that these polarons become increas- 
ingly correlated with each other, producing short range 
charge ordering on lowering the temperature towards Tc- 
Both the quasistatic polaronic scattering and the charge 
correlation peaks start to collapse just above Tc, and dis- 
appear in the ferromagnetic, metallic state (see Fig. H). 
It has been argued that these polarons explain the low 
hole mobility in the paramagnetic state, which cannot be 
due to double exchange alone [Q. 

The present observations are strong evidence that the 
phase transition in Lai.2Sri.gMn207 is driven by the 
growth of magnetic correlations, which are only weakly 
coupled to the polarons above Tc- The polarons are 
likely to induce some exchange disorder, but it is evi- 
dently not sufficient to disrupt the development of critical 
magnetic fluctuations. This conclusion is not consistent 
with a magnetic polaron model, in which the spin corre- 
lations are strongly coupled to the charge degrees of free- 
dom. In such models, the ferromagnetic phase transition 
is induced by a transition from small to large polarons 
1^ , but it is unlikely that such a transition would mimic 
the scaling behavior observed here. This docs not mean 
that such models cannot be relevant to other CMR com- 
pounds, particularly those with much stronger electron- 
phonon coupling, but they are not essential to a descrip- 
tion of the CMR process. 

Although the spin-charge coupling is weak above Tc, 
the magnetic correlations ultimately drive the metal- 
insulator transition; once the spin correlations extend 
over a large enough region, the double exchange inter- 
action can overcome the mechanism responsible for lo- 
calizing the charges, inducing the polaron collapse [^. 
This may affect the critical scaling in the 3D regime and 
could be responsible for the unusually low value of (3 (in 
quasi-2D XY systems, f3 is predicted to be approximately 
0.23 0). 

Our results demonstrate that the critical properties 
of the bilayer CMR manganite Lai.2Sri.8Mn207, both 
above and below Tc , are in quantitative agreement with 
an effective finite-size 2D XY model. The magnetic and 
charge degrees of freedom are therefore only weakly cou- 
pled except close to Tc, where the growth of magnetic 
order delocalizes the polaronic charges. It is possible that 
some form of charge ordering would occur at lower tem- 
perature if it were not preempted by the magnetic order- 
ing. From these observations, we conclude that magnetic 
polaron models are not appropriate to the present bilayer 
compound and are therefore not universal to CMR. 
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